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Abstract 

We study the normal state conductance and the Josephson current in a 
superconductor-2DEG-superconductor structure where the size/shape of 
the 2DEG-region can be modified by an additional side-gate electrode. 
The considered transport properties follow from the retarded Green func- 
tion which we compute by employing a tight-binding-like representation 
of the Hamiltonian in the 2DEG region. Our model studies offer a quali- 
tative demonstration of the recently observed effects caused by side-gate 
modulation. 
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1. Introduction 

Devices with semiconducting materials as barriers in Josephson weak links of- 
fer the possibility to modulate the junction properties in order to investigate 
both fundamental transport mechanisms as well as for use in applications of 
superconducting devices. Several methods for tunable Josephson junctions have 
been studied both theoretically and experimentally. Tunable effects have been 
observed in e.g. Josephson-field-effect-transistors [Q, non-equilibrium junctions 
0, 3, optically modulated weak links 0, and recently in side-gate modulated 
junctions p. 

Following Ref. we consider superconductor-2DEG-superconductor struc- 
tures where the indirect ballistic transport between the non-opposite supercon- 
ducting contacts can be controlled by a voltage applied to a side-gate (see Fig. 



1 
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0). The geometry has many similarities with the T-stub wave guide geometry 
1^, 1^ but here the two interfaces to the superconductors are non-opposite. The 
aim is to demonstrate the principle behind side-gate modulation of the critical 
current and also to investigate the local density-of-states in the 2DEG in relation 
to the discussion of "non-local" modes ^. 

In this work we focus on the normal state conductance and the Josephson 
current in the case where electron and hole propagation in the normal region 
is phase- coherent. For convenience we consider the zero temperature limit. The 



normal state conductance is then given by the Landauer formula ||10| , |Tl|, [12 
which at T = K has the particularly simple form 

GN = ^Trttt = ^^T., (1) 

n=l 

where 2e^/h is the quantum unit of conductance and T„ (n = 1, 2, . . . A^) are the 
transmission eigenvalues of tt^ which is an x matrix, being the number of 



propagating modes at the Fermi level. Beenakker and van Houten ||T3[ considered 
the Josephson current through a superconducting quantum point contact and 
found that the critical current was quantized in units of eA/h, A being the 



energy gap of the superconductor. Subsequently it was shown by Beenakker |T4 
that just like the normal state conductance, also the Josephson current between 
two superconducting leads can also be expressed in terms of of the normal 
region which couples the two superconductors. At T = K the Josephson current 
can be written as [ p!^ 



= sin 0Tr tt^ \l - tt^ sin^ (0/2)1 ^ 
2n -' 

= |^sin0^T„[l-T„sin2(0/2)]-'/', (2) 

n=l 

where is the phase difference between the two superconductors. Eq. is 
valid in the short-j unction limit where the dimension of the normal region is 
short compared to the superconducting coherence length ^ which in the ballistic 
regime is ^ = Hvp/ttA, vp being the Fermi velocity. 

We note that in Eqs. (|I]) and (0) all transverse degrees-of-freedom are incor- 
porated in the transmission amplitude matrix t and that the angle- dependence 
of the Andreev scattering discussed in Ref. UTS follows directly from the an- 



gle (mode) dependence of t. Scattering due to non-matching Fermi velocities or 
Fermi momenta can be included in several ways: i) hj calculating the transmis- 
sion from the composite scattering matrix S = Si^^ ® 5'2deg 8) 5*12 where the two 
interfaces (Ii and I2) are described by the scattering matrices 5*1^ and 5*12 as in 
Ref. |]T6[, ii) by adding an effective potential to the 2DEG at the interfaces, or 
Hi) by explicitly taking the different Fermi velocities and Fermi momenta into 
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Figure 1: Plot of Gn, ^c, and /c-Rn as a function of transmission probability T = |ip for a 
system with single-mode leads. 

account. In this work we will however for simplicity assume matching Fermi 
properties of the superconductor and the 2DEG. 

To illustrate how the transport properties depend on the transmission let us 
consider the case of single-mode leads. Fig. ^ shows how the normal state con- 
ductance Gn, the critical current Ic = max/j(0), and the product /c-Rn = -^c/G^n 
depend on the transmission probability T = Compared to the off- resonance 
regime with T ~ (the Ambegaokar-Baratoff regime [1^) the /c-Rn product 
is enhanced by a factor-of-two when the 2DEG region is tuned to resonance 
(T ~ 1). This factor-of-two enhancement can be considered as a signature of 
Andreev mediated transport at resonance whereas the transport is tunneling- 
like at off-resonance. 

Equations ([l|) and (||) form the basis for our calculations of the transport prop- 
erties and since they only depend on the normal state transmission properties 
of the 2DEG region we can employ standard methods for quantum transport in 
semiconductor structures. The transmission can conveniently be calculated from 
the retarded Green function ^J(ri,r2) ||12[ and also the local density-of-states 
follows from ^J(ri,r2). We calculate ^J(ri,r2) by applying a finite differences 
method to the Hamiltonian of the 2DEG region. 

The paper is organized as follows: Section II introduces the finite differences 
method and formulates Eqs. (|l|) and (0) in terms of the retarded Green function, 
section III presents our geometry and model, section IV contains our results, and 
in section V discussions and conclusions are given. 




2. Finite differences method 

We discretize the continuous Hamiltonian by writing the Laplacian as finite 



differences 18 . This gives rise to a tight-binding-like representation 



{i/eh, = zX7 + [/(r,) t=j, (3) 
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= -7 ^,JNN, 
= otherwise, 

where 2; = 4 is the number of nearest neighbors (NN) and 7 = fi^ /2ma? corre- 
sponds to a hopping matrix element, a being the lattice spacing. Comparing the 
energy dispersion to the parabolic dispersion of the continuous problem shows 
that the finite differences method is a very accurate description for energies e < 7 
T8| , [T9[| . Making the grid finer (a smaller) improves the accuracy and/or allows 



for a treatment of higher energies. 

With the above definition of the Hamiltonian we label the lattice points by a 
number from 1 to M, M being the number of lattice points. The retarded Green 



function can then be written as an M x M matrix |T^, |T9 

g\e)=[el~H,-T.\e)Y\ (4) 

where the element {Q^{£)}ij corresponds to ^J(rj,rj). Here, He is the tight- 
binding Hamiltonian of the conductor (describing the conductor as a closed 
system) and = T,\ + is the self-energy describing the coupling to leads 1 
and 2. For the self-energies the only non-zero elements are those where the leads 
couple to the 2DEG (scattering region). For large M it can be useful to employ a 
recursive method for calculating Q^{e) but even for M ~ 1000 the direct matrix 
inversion in Eq. (^ can be done without too much numerical effort. 

The transmission amplitude matrix t can now be calculated from the Fisher- 
Lee relation [0 which gives 0, |TH|, 0, 

t{s) = [T,{e)f'G\e)[T,{e)f\ (5) 
where T j{e) = i (^^(e) - m{e)]^). Substituting into Eqs. d) and (D gives 

= —TTT,{e)g^-{e)T,{eWie)]\ (6) 

and 



IM) = |^sin0Trri(e)^>^(.)r2(5)r(5)]t (7) 



-1/2 



where we have used the cyclic invariance of the trace. Once the retarded Green 
has been obtained (by a single matrix inversion) these relations directly provide 
the essential transport properties. The local density-of-states (in the normal 
state) can also be calculated directly from the retarded Green function (see e.g. 

m) 

P(r„^) = -^Im{^'^(e)},,, (8) 
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Figure 2: Indirect S-2DEG-S structure with a side-gate (left) and the corresponding lattice 
model (right). 

and this function is useful in obtaining insight into the spatial variations of the 
states and revealing the nature of the conduction through the sample. 

3. Geometry and model 

We consider the geometry shown in Fig. 0. The 2DEG is confined also in the 
direction parallel to the side- gate so that the system acts as a cavity or quantum 
dot coupled to two superconductors and an additional side-gate. We note that our 
geometry is slightly different from that of Ref. where there is no confinement 
of the 2DEG in the direction parallel to the side-gate. This means that here the 
side-gate is used in tuning the cavity to resonance whereas in Ref. it rather 
acts as a 'classical' mirror which can be used in focusing the incident wave from 
one lead onto the other lead. 

We treat the transmission problem fully quantum mechanically by means of 
the presented finite differences method and the lattice version of the sample 
shown in Fig. ^ In this case a = W/ (Ml + 1) with Ml = 4 sites in the transverse 
direction of the leads of width W. In the continuous description of the leads 
= Intlk-FW/ir) gives the number of propagating modes at the Fermi level, 
Int(x) being the integer part of x. The dimensions are 

and W" = (13 + l/2)V2a ~ 3.8 x W. We choose W and the Fermi level such that 
N = 1. The lattice model then gives a reasonable description of the continuous 
problem - the threshold energy Ei of the first mode deviates only by ~ 3% from 
the continuous result. 

The side-probe is assumed to act as a gate and only affect the potential U of 
the 2DEG through an electro-static coupling. Here, we note that fluctuations in 
the gate potential may lead to dephasing of the electron and hole propagation in 
the 2DEG (see e.g. Ref. [^) and the same can also be the case if the probe acts 




Figure 3: The transmission of the lowest mode T = |ip as a function of the energy e (nor- 
malized by the threshold energy E\) for gate voltages = j x 8V^ with = 0, 1, 2, . . . 10. 



as a voltage probe (see e.g. Ref. [|^). Here, we use the model in Ref. [Eq. 
(19)] to account for the potential modification in the 2DEG due to the side-gate. 
The distance C from the side-gate to the 2DEG edge position is then given by 

eVg eVg , . 

^ " 477% |e| " 27rfc2|e|' ^ ^ 

where is the side-gate voltage, e = e^Co is the dielectric constant of the semi- 
conductor (e,- ~ 14 for GaAs), and = kp/27T the 2DEG density. With the 
lattice shown in Fig. ^ the distance C can be changed in steps of 6C = a/V2 
corresponding to 5Vg = ^/2nkpa\e\/e. 



4. Results 

In Fig. ^ we show the transmission as a function of the energy for different 
values of the gate voltage. Since the interfaces between the leads and the 2DEG 
are fully transparent the transmission shows a broad resonance behavior. By 
changing the gate voltage the cavity is 'squeezed' and as expected the resonance 
shifts towards higher energies. For a given Fermi level the side-gate can thus also 
be used in tuning the transmission to resonance. 

In Fig. ^ we show the normal state conductance G^, the critical current Jc, 
and the /c-Rn product as a function of gate voltage for an energy e = 2.5 x Ei. 
By increasing the gate voltage the cavity is tuned to resonance at ~ 6 x 6Vg 
where the normal state conductance equals the quantum unit of conductance 
2e^ /h and the critical current equals the quantum unit of critical current eA/^. 
The /c-Rn product changes by a factor-of-two from the Ambegaokar-Baratoff 
value A/87re to its quantum unit A/47re at resonance. 
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Figure 4: The normal state conductance Gn, the critical current Ic, and the IcRn product as 
a function of gate voltage Vg for an energy e = 2.5 x E'l. The inset shows the local density of 
states p{r, e) (in the normal state) at gate voltage Ig = 6 x SVg corresponding to the resonance 
condition. Dark regions indicate a high density of states {e.g. near the leads) and bright regions 
a lower {e.g. near the side-gate). 



In the inset of Fig. § we show the local density of states p(r, e) (in the normal 
state) at the energy e = 2.5xEi for a gate voltage corresponding to the resonance 
condition. The transparent interfaces give rise to a local density of states forming 
a relatively smooth and connected 'path' between the two leads. The plot also 
clearly shows the wave nature of the electron and hole propagation. Thus a more 
simple semi-classical trajectory model with point-like particles and a side-gate 
acting as a classical mirror would be inadequate for the present situation. This 
is often the case and as studied in e.g. Ref. complicated multiple scattering 
processes can in a full quantum mechanical treatment give rise to effects not 
found in a semi-classical study. 



5. Discussion and conclusion 

Hybrid semiconductor-superconductor structures offer interesting possibilities 
for investigating fundamental transport phenomena as well as for potential appli- 
cations. We have studied the possibility of side-gate modulating the critical cur- 
rent in a mesoscopic superconductor-2DEG-superconductor Josephson junction. 
Side-gate modulation offers a new alternative to devices based on field-effects, 
non-equilibrium effects, and optical effects. In the case of side-gate modulation 
the side-gate is used in tuning the transmission of the 2DEG to resonance. 
Our calculations are based on a numerical treatment of the transmission prop- 



erties of the 2DEG region. Using the Landauer formula [T^, |rT| |T2[ and a similar 



formula for the Josephson current we have then calculated the essential 
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transport properties: the normal state conductance Gn, the critical current /c, 
and the /c-Rn product. The conductance in the superconducting state can unfor- 
tunately not be expressed directly in terms of the transmission eigenvalues and 
is complicated due to the presence of multiple Andreev reflections (see however 
Ref. for a possible solution). The considered geometry is comparable to the 
situation in recent experiments even though some simplifications have been 
made. A quantitative comparison would require i) that the different Fermi prop- 



erties of the superconductor and the 2DEG are taken into account [|T5], ii) more 
lattice points (Ml 2> N) in order to account correctly for the case of multi- 
mode leads (A^ ~ 16 in Ref. 0), and Hi) that the width of the 2DEG should be 
much larger than the separation W of the two leads. However, our model studies 
demonstrate how the side-gate modulation can be used in controlling both the 
normal state conductance and the critical current. 

Our results are in qualitative agreement with the recent experimental findings 
in and confirm the possibility of modulating the transport properties by a side- 
gate. Studies of the local density of states confirm the existence of "non-local" 
modes and indicate that more simple trajectory models can not fully account for 
the detailed electron and hole propagation - rather a full quantum mechanical 
treatment is necessary. 
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